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Abstract. Globally mapped terrestrial chlorophyll fluores- 
cence retrievals are of high interest because they can provide 
information on the functional status of vegetation includ- 
ing light-use efficiency and global primary productivity that 
can be used for global carbon cycle modeling and agricul- 
tural applications. Previous satellite retrievals of fluorescence 
have relied solely upon the filling-in of solar Fraunhofer lines 
that are not significantly affected by atmospheric absorption. 
Although these measurements provide near-global coverage 
on a monthly basis, they suffer from relatively low preci- 
sion and sparse spatial sampling. Here, we describe a new 
methodology to retrieve global far-red fluorescence infor- 
mation; we use hyperspectral data with a simplified radia- 
tive transfer model to disentangle the spectral signatures of 
three basic components: atmospheric absorption, surface re- 
flectance, and fluorescence radiance. An empirically based 
principal component analysis approach is employed, primar- 
ily using cloudy data over ocean, to model and solve for 
the atmospheric absorption. Through detailed simulations, 
we demonstrate the feasibility of the approach and show 
that moderate-spectral-resolution measurements with a rela- 
tively high signal-to-noise ratio can be used to retrieve far-red 
fluorescence information with good precision and accuracy. 
The method is then applied to data from the Global Ozone 
Monitoring Instrument 2 (GOME-2). The GOME-2 fluores- 


cence retrievals display similar spatial structure as compared 
with those from a simpler technique applied to the Green- 
house gases Observing SATellite (GOSAT). GOME-2 en- 
ables global mapping of far-red fluorescence with higher pre- 
cision over smaller spatial and temporal scales than is possi- 
ble with GOSAT. Near-global coverage is provided within a 
few days. We are able to show clearly for the first time phys- 
ically plausible variations in fluorescence over the course of 
a single month at a spatial resolution of 0.5° x 0.5°. We also 
show some significant differences between fluorescence and 
coincident normalized difference vegetation indices (NDVI) 
retrievals. 


1 Introduction 

Vegetation releases unused absorbed photosynthetically ac- 
tive radiation primarily as heat with a small amount re- 

emitted as fluorescence. Measurements of terrestrial chloro- 
phyll fluorescence are directly related to photosynthetic 

function, and are potentially useful for forest and agri- 
cultural applications as well as assessment of the terres- 
trial carbon budget by providing more accurate estimates of 
gross primary productivity (GPP) (e.g., Lichtenthaler, 1987; 

Saito et al., 1998; Corp et al., 2003, 2006; Campbell et al., 
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2008; Damm et al., 2010; Joiner et al., 2011, 2012; Franken- 
berg et al., 2011b; Guanter et al., 2012). Studies show that in 
high light conditions, such as in the late morning and early 
afternoon, when many satellite measurements are made and 
when plants are under stress, fluorescence is correlated to 
the amount of absorbed photosynthetically active radiation 
(APAR) and the efliciency of the plants to utilize this light 
to drive photosynthesis (light-use efliciency or LUE) (e.g., 
Flexas et al., 2002; Louis et al., 2005; Meroni et al., 2008; 
Amoros-Lopez et al., 2008; van der Tol et al., 2009; Zarco- 
Tejada et al., 2009, 2013; Daumard et al., 2010). 

Fluorescence information is also complementary to 
reflectance-based spectral vegetation indices (Meroni and 
Colombo, 2006; Middleton et al., 2008, 2009; Rascher et al., 
2009; Meroni et al., 2008; Daumard et al., 2010; Guanter 
et al., 2007, 2012; Zarco-Tejada et al., 2009, 2013; Joiner 
et al., 2011, 2012; Frankenberg et al., 2011b). For example, 
popular greenness-based indices such as the normalized dif- 
ference and enhanced vegetation indices (NDVI and EVI, re- 
spectively) are linked to chlorophyll content and are related 
to potential photosynthesis, whereas fluorescence should be 
an indicator of actual photosynthesis. The photochemical re- 
flectance index (PRI) is sensitive to the de-epoxidation state 
of pigments within xanthophyll cycle, a protection mech- 
anism that evolved in parallel to fluorescence to dissipate 
excess energy (Gamon et al., 1992). 

Knowledge of global chlorophyll fluorescence emissions 
is also important for retrievals of trace-gas concentrations 
including CO 2 that require very high accuracy and preci- 
sion. The emission occurs within the O 2 A-band that is 
used to estimate photon path lengths for these measure- 
ments. If not properly accounted for, fluorescence emis- 
sion may produce significant errors in these retrievals 
(Frankenberg et al., 2012). The effect of fluorescence on an 
aerosol plume height retrieval from the O 2 A-band has also 
been investigated using a linear error analysis with simulated 
data (Sanders and de Haan, 2013). 

One means of measuring the small fluorescence signal 
from passive remote sensing instrumentation is to make use 
of dark features in the Earth’s reflected spectrum, either from 
telluric absorption or deep solar Fraunhofer lines. For ex- 
ample, ground-, aircraft-, and space-based approaches have 
utilized filling-in of the dark and spectrally wide O 2 A- 
band (~760nm) and O 2 B-band (~690nm) atmospheric 
absorption features to detect the weak fluorescence signal 
(see, e.g., Guanter et al., 2007; Meroni et al., 2009, 2010). 
The spectral location of these oxygen absorption features as 
well as other absorption bands and solar Fraunhofer lines 
are shown in Fig. 1 along with the broadband red and far- 
red fluorescence emission features that peak near 685 and 
740 nm, respectively. 

Deep solar Fraunhofer lines have also been used to de- 
tect fluorescence from vegetation following the early work 
of, e.g., Plascyk and Gabriel (1975). Joiner et al. (2011, 
2012), Frankenberg et al. (2011b), and Guanter et al. (2012) 



Wavelength (nm) 

Fig. 1 . Simulated solar- induced terrestrial fluorescence, typical sim- 
ulated atmospheric transmittance and reflectance, and solar irradi- 
ance as a function of wavelength computed for an instrument with 
FWHM = 0.3 nm. The fluorescence shows red and far-red chloro- 
phyll emission features with peaks near 685 and 740 nm, respec- 
tively. Oxygen A and B absorption bands are located near 687 and 
760 nm, respectively, while water vapor absorption is shown over 
a broad spectral range between about 690 and 740 nm. The solar ir- 
radiance shows weak solar Fraunhofer line structure at this spectral 
resolution. 

used near-infrared (NIR) solar Fraunhofer lines, which are 
filled in by vegetation fluorescence, to globally map terres- 
trial fluorescence with the high-spectral-resolution interfer- 
ometer aboard the Japanese Greenhouse gases Observing 
SATellite (GOSAT). Joiner et al. (2012) strongly suggested 
that fluorescence may also be measurable from space with 
lower-spectral- resolution instrumentation as compared with 
the GOSAT interferometer or similar instruments. They fo- 
cused on filling-in of the 866 nm Ca II solar Fraunhofer line 
as measured with the SCanning Imaging Absorption spec- 
troMeter for Atmospheric CHartographY (SCIAMACHY) 
satellite instrument. This filling-in appears to be produced by 
fluorescence from chlorophyll as supported by Gamon and 
Berry (2012). The SCIAMACHY spectral resolution at this 
wavelength is about 0.5 nm. 

While the current satellite results show promise for use in 
estimation of GPP, the GOSAT measurements have fairly low 
spatial sampling and relatively low single-observation preci- 
sion (Joiner et al., 2011; Frankenberg et al., 201 la; Guanter 
et al., 2012). The SCIAMACHY results have higher sam- 
pling frequency, but the very low signal levels spectrally far 
from the far-red fluorescence peak also result in low pre- 
cision for single observations (Joiner et al., 2012). To pro- 
duce global maps with high enough fidelity for comparisons 
with other measurements and models, GOSAT and SCIA- 
MACHY fluorescence retrievals must be averaged spatially 
and/or temporally. In doing so for GOSAT, there can be 
a substantial sampling or representativeness error introduced 
by the averaging of sparse observations within a relatively 
large grid box. 
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Other approaches for satellite fluorescence retrievals have 
aimed at utilization of the strong atmospheric oxygen bands 
(A and B bands) that absorb at wavelengths where chloro- 
phyll fluorescence is emitted. For example, approaches to 
separate fluorescence features from those of reflectance for 
space-based measurements have been developed in which at- 
mospheric absorption is assumed to be perfectly modeled 
(e.g., Mazzoni et al., 2008, 2010, 2012). More complex al- 
gorithms have been proposed and tested on simulated data in 
which the parameters affecting O 2 absorption are retrieved 
and accounted for using a radiative transfer model (Guan- 
ter et al., 2010). Low- spectral-resolution O 2 A-band satellite 
measurements from the MEdium Resolution Imaging Spec- 
trometer (MERIS) aboard Envisat have also been used to re- 
trieve information about fluorescence (Guanter et al., 2007). 
Thus far, these satellite measurements are spatially and tem- 
porally limited and require an on-ground nearby nonfluoresc- 
ing reference target for normalization. 

Here, we develop new methodology to retrieve the far- 
red chlorophyll fluorescence using space-based hyperspec- 
tral measurements in and around the O 2 A-band. Instead of 
exclusively using the filling-in of solar Fraunhofer lines as 
in the previous works with GO SAT and SCI AM ACHY, we 
demonstrate that fluorescence can be retrieved by exploit- 
ing the different spectral structure produced by the far-red 
chlorophyll fluorescence feature (including both solar and 
telluric line filling), atmospheric absorption, and surface re- 
flectance. 

Our methodology is similar to approaches developed for 
ground-based instrumentation (Guanter et al., 2013) in that 
radiative transfer in atmospheric absorption bands is ap- 
proximated statistically using a principal component anal- 
ysis (PC A) (or singular value decomposition, SVD). SVD 
approaches have also been applied to satellite fluorescence 
retrievals using wavelengths not affected by atmospheric ab- 
sorption (Guanter et al., 2012). In this work, we expand the 
use of PCA to include the geometry of a space-based in- 
strument for wavelengths where significant atmospheric ab- 
sorption takes place. This scenario is more complex than for 
ground-based measurements. Fluorescence emission can be 
significantly absorbed in the atmosphere. Because this ab- 
sorption is different in magnitude from that of reflected sun- 
light, the scenario for a satellite retrieval is more difficult as 
compared with that for a ground-based instrument. While our 
approach does not require a nearby nonfluorescing target as 
in other techniques, a representative sample of observations 
over nonfluorescing scenes is needed in order to generate a 
comprehensive set of PCs. For this purpose, we use cloudy 
observations over ocean covering a large range of latitudes 
on a daily basis. 

We conduct simulations to demonstrate the applicability of 
our approach to current and proposed satellite instruments. 
We then apply our technique to data from the Global Ozone 
Monitoring Instrument 2 (GOME-2). The primary function 
of GOME-2 is to make measurements of atmospheric trace 
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gases. While not optimal for fluorescence measurements ow- 
ing to its relatively large ground footprint and moderate spec- 
tral resolution, its high sampling and signal-to-noise ratio 
enable state-of-the-art fluorescence retrievals in the far-red 
chlorophyll emission feature. Near-global coverage is pro- 
vided within a few days from GOME-2 measurements with 
a high signal-to-noise ratio. Accurate and frequent measure- 
ments from GOME-2 will lead to a fluorescence data set with 
unprecedented temporal and spatial resolution; this in turn 
should enable detailed studies including more direct compar- 
isons with flux tower measurements. 

2 GOME-2 satellite data 

In this work, we use data from GOME-2. GOME-2 is an 
operational nadir-viewing UV/visible cross-track scanning 
spectrometer (Munro et al., 2006). It flies as part of the Eu- 
ropean Meteorological Satellite (EUMETSAT) Polar System 
(EPS) MetOp mission series. GOME-2 measures the Earth’s 
backscattered radiance and the extraterrestrial solar irradi- 
ance at wavelengths between 240 and 790 nm in four detec- 
tor channels. Here, we use level IB data from revision R2 in 
channel 4 that cover wavelengths 590-790 nm with a spectral 
resolution of approximately 0.5 nm (Callies et al., 2000) and 
a relatively high SNR (> 1000). 

The nadir Earth footprint size is 40 km x 80 km, and the 
nominal swath width is 1920 km. A single GOME-2 in- 
strument provides global coverage of the Earth’s surface in 
about 1.5 days. The first flight of GOME-2 is on MetOp- A, 
launched 19 October 2006 into a polar orbit with an equa- 
tor crossing time of 09:30 LT. The second flight, launched 
17 September 2012 on MetOp-B, is also in a morning orbit 
but 180° out of phase with the first flight model. As such, 
one or the other instrument is always making measurements 
on the sunlit part of the Earth, and near-global daily coverage 
is achievable. 

3 Simulated radiances and irradiances 

Because our approach relies on an empirical rather than 
a physical approach for deriving atmospheric absorption, it 
is difficult to quantify forward model errors. To accurately 
quantify retrieval errors, we conduct detailed simulations 
using combined atmospheric and vegetation models over 
a wide range of conditions. We also use simulated data to 
assess the impact of instrument specifications including the 
signal-to-noise ratio (SNR) and spectral resolution on fluo- 
rescence retrievals. Finally, we test different retrieval scenar- 
ios, such as various spectral fitting windows and numbers of 
retrieved parameters, using the simulated radiances. 

We simulate top-of-the-atmosphere (TOA) sun- 
normalized radiances using the matrix operator model 
(MOMO) radiative transfer model (Fell and Fischer, 2001; 
Preusker and Lindstrot, 2009). The radiance calculations 
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Fig. 2. Simulated solar spectra based on Chance and Kuracz (2010) 
for different instrumentation showing Fraunhofer line structure 
(FWHM = 0.5 nm is for a GOME-2-like instrument). 
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Fig. 3. Canopy-level spectral fluorescence as specified for the range 
of conditions in the simulated data set. 


utilize absorption line strengths and widths from the high- 
resolution atmospheric radiance and transmittance model 
code (HITRAN) 2008 data set (Rothman et al., 2009). The 
radiances are computed monochromatically and are sampled 
at 0.005 nm. They are then multiplied by a solar spectrum 
sampled in the same way and finally convolved with various 
instrument line shape functions. 

We use solar data originally sampled at 0.001 nm from 
kurucz.harvard.edu/sun/irradiance2005/irradthu.dat, similar 
to Chance and Kurucz (2010) but more highly sampled. Fig- 
ure 2 shows simulated solar spectra generated for different 
instrument specifications including a spectral resolution sim- 
ilar to GOME-2 (i.e., full width at half maximum (FWHM) 
of 0.5 nm sampled at 0.2 nm) and a smaller FWHM of 0.3 nm 
sampled at 0. 1 nm. Significant solar Fraunhofer line structure 
can be seen throughout the spectrum with deeper structures 
at the higher spectral resolution. 

Note that we do not simulate the effects of rotational 
Raman scattering (RRS) or O 2 A-band day glow emissions. 
RRS effects are generally small, though not negligible, at the 
wavelengths of interest (Vasilkov et al., 2013). The effects of 
O 2 A-band dayglow emissions in the upper atmosphere are 
expected to be small (Guanter et al., 2010). Directional ef- 
fects of the vegetation reflectance and fluorescence are also 
not simulated. 

Radiances are computed for two view zenith angles 
(VZA = 0 and 16°), four solar zenith angles (SZA=15°, 
30°, 45°, 70°), two atmospheric temperature profiles (mid- 
dle latitude summer and winter), four surface pressures (955, 
980, 1005, and 1030hPa), four values of total column wa- 
ter vapor (0.5, 1.5, 2.5, and 4.0 cm), five aerosol optical 
thicknesses at 550 nm (0.05, 0.12, 0.2, 0.3, and 0.4), and 
three aerosol plume heights (700-900, 600-800, and 500- 
700 hPa) with a continental aerosol model. Two separate data 
sets are created, one without fluorescence intended for prin- 


cipal component analyses (henceforth referred to as “train- 
ing”), and one containing fluorescence intended to examine 
retrieval performance with the simulated radiances (referred 
to as “testing”). The training data set uses a spectral library 
of 10 different surface reflectance spectra corresponding to 
soil and snow. This provides a total of 38 400 samples in the 
training data set. 

The testing data set contains surface reflectance and flu- 
orescence emission spectra generated with the FluorSAIL 
and FluorMODleaf codes (Jacquemoud et al., 2009; Pe- 
dros et al., 2010; Miller et al., 2005). A given input leaf- 
level fluorescence spectrum has been scaled by multiplica- 
tive factors. Different values of chlorophyll content (from 
5 to 40 pg cm -2 ) and leaf area index (from 0.5 to 4) have 
been used for the propagation of the resulting leaf-level flu- 
orescence spectra to the canopy level. There are 60 distinct 
top-of-canopy fluorescence spectra from these combinations 
as shown in Fig. 3. Top-of-canopy reflectance spectra are 
consistently generated by the leaf and canopy codes for the 
same combinations of chlorophyll content and leaf area in- 
dex. Other parameters in the models were set to default val- 
ues as in Guanter et al. (2010). This gives a total of 230400 
different samples in the testing data set. 


4 Retrieval methodology 
4.1 General approach 


The total reflectance pt 0 t, as a function of wavelength k, mea- 
sured by a satellite instrument in the NIR spectral region, can 
be approximated using a Lambertian model with emission 
from fluorescence, i.e., 


PtotG) — poG) + 


p s (k)T(k)T(k) 

l-p s (A)p(A) 


7 tF s (X)T(X) 

[l-p s (A)p(A)]F(A)cos(P 0 ) a } 
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where p s (A) is the surface reflectance, po is the reflectance 
contribution in the absence of surface effects, T is the total 
irradiance transmittance (including direct and diffuse com- 
ponents), T is the spherical transmittance from the surface to 
TOA, p is the spherical reflectance of the atmosphere back to 
the surface, 0o is the solar zenith angle (SZA), E(X) is the ob- 
served extraterrestrial solar irradiance, and F s is the radiance 
emission from fluorescence at the surface. 

The basic idea behind our approach is to separate the spec- 
tral features related to three basic components: atmospheric 
absorption (T and T), surface reflectivity (p s ), and fluo- 
rescence radiance ( F s ). Assuming that the effects of atmo- 
spheric scattering are small (i.e., po — 0 and p s p 1), we 
may rewrite Eq. (1) as 

AotW = p^TixyTiX) + (2) 

£(A)cos(0o) 


smooth, and model it as a low-order polynomial in A. Al- 
ternative parameterizations for fluorescence and reflectance 
have been explored (e.g., Mazzoni et al., 2010, 2012). 
Previous works suggest that small errors in the prescribed 
shape of the fluorescence emission have little impact on the 
estimated peak fluorescence value (Daumard et al., 2010; 
Fournier et al., 2012; Guanter et al., 2013). We estimate the 
spectral structure of A 2 (or T 2 ) using principal components 
(PCs) as described below. 

In principle, our approach may be applied to the entire flu- 
orescence emission band shown in Fig. 1 containing both the 
red and far-red features. Alternatively, different fitting win- 
dows could be used to estimate fluorescence within smaller 
wavelength ranges. As a starting point to demonstrate our 
approach, we focus on retrievals of the far-red fluorescence. 

4.2 Generation of atmospheric PCs 


where T (A) and T (X) now include atmospheric molecu- 
lar absorption only. We now combine T and T into a sin- 
gle parameter, 72(A) = T(A)T(A ), that represents the sun- 
to-satellite (two-way) atmospheric transmittance. Using 

T 2 {X) = exp [—*42 (A)] 

= exp[-*4 v (A){sec(0) + sec(0 o )}], (3) 


where 0 is the satellite view zenith angle (VZA), and *42 and 
A v represent the two-way and vertical absorptances, respec- 
tively, and the upwards absorptance, *4, as 


A(X) = A v (A) sec (0) 

sec ( 0 ) 




sec (0) + sec (0 q) ’ 


then 


(4) 


T (A) = exp ( In [72 (A)] 


sec (0) 


sec (0) + sec (0 q) 


(5) 


Note that the above equations are strictly valid only for 
monochromatic radiation. Also note that po, p, T, and T, 
in the presence of atmospheric scattering and in the absence 
of atmospheric molecular absorption, are a spectrally smooth 
function of wavelength. Therefore, when atmospheric scat- 
tering is present, p s (A) and F S (A) in Eq. (2) can be thought 
of as effective TOA spectral components of surface re- 
flectance and fluorescence that have been modified by spec- 
trally smooth atmospheric scattering; the spectral structure of 
po can be incorporated into the components of the first term 
of Eq. (2). 

To solve for p s , F s , and T 2 , we assume that 
each has a distinct spectral structure. We repre- 
sent the fluorescence far-red emission, F S (A), as 

a Gaussian function of A centered at 736.8 nm with 
a = 21.2 nm, similar to Subhash and Mohanan (1997) and 
Zarco-Tejada et al. (2000). We further assume that p s (A), 
within our limited spectral fitting window, is spectrally 


Radiative transfer in the O 2 A-band is complex because the 
shape of the absorption features depends nonlinearly upon 
the surface pressure and albedo as well as the vertical struc- 
ture of atmospheric temperature and cloud/aerosol particles 
(e.g., Preusker and Lindstrot, 2009). Absorption by water va- 
por in the 710-745 nm spectral region also depends upon 
these parameters and is similarly complex; though less af- 
fected by saturated lines, the profile is of course more vari- 
able. Retrieval algorithms have been developed to retrieve 
some or all of these parameters using O 2 A-band radiances 
(e.g., Guanter et al., 2007; O’Dell et al., 2012). Other factors 
such as filling-in from rotational Raman scattering are typi- 
cally not accounted for owing to a high degree of complexity 
resulting from interactions with Mie scattering and depen- 
dences on other parameters (Vasilkov et al., 2013). Neglect of 
these effects may produce errors in such algorithms. Instru- 
mental effects such as nonlinearity also challenge physically 
based approaches used to estimate O 2 A-band absorption. 

Instead of using radiative transfer calculations, we have 
developed an empirically based alternative for estimation of 
*42 or 75(A); we represent *42 (A) as a linear combination of 
principal components (PCs) fa (A) that can be estimated us- 
ing simulated or real satellite data, i.e., 

n 

A 2 (X) = J2 a '^ (X) ’ (6) 

i = 1 

where ai are the coefficients of the PCs. Instead of using 
laboratory-measured absorption cross sections, as is typical 
in the differential optical absorption spectroscopy (DOAS) 
approach, we use atmospheric spectra, simulated or mea- 
sured, to derive the spectral components of absorption. 

As in the DOAS approach, our method implicitly assumes 
that the Beer-Lambert law of weak linear absorption applies, 
although the PCs may be able to incorporate some features 
of nonlinear absorption. Because this law does not strictly 
apply to the O 2 A-band, where individual lines may become 
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Fig. 4. Transmittances derived from the training data set with FWHM = 0.5 nm. Colors in the left panel correspond to values of total column 
water vapor (light blue, dark blue, green, and red correspond to 4.0, 2.5, 1.5, and 0.5 cm, respectively) for various combinations of solar and 
viewing zenith angles; colors on right-hand panel showing O 2 A-band absorption correspond to different combinations of solar and viewing 
zenith angles. 


0.15F 
0.10 
0.05 
0.00 


PC % of total: 99.8392(99.6467); cumulative % of total: 99.8392(99.6467) 



715 720 725 730 735 

Wavelength (nm) 


740 


PC % of total: 0.1301 ( 0.3244); cumulative % of total: 99.9693(99.971 1) 



745 


o 

Q_ 


0.2 - 
0.1 
0.0 
- 0.1 


0.2 F - 
0.1 - 
0.0 H 

ML 


715 720 725 730 735 

Wavelength (nm) 


740 


PC % of total: 0.0224( 0.0198); cumulative % of total: 99.9917(99.9910) 


PC % of total: 0.0041 ( 0.0064); cumulative % of total: 99.9958(99.9974) 


745 



o 

Q_ 



715 720 725 730 735 

Wavelength (nm) 


740 


745 
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Fig. 6. Similar to Fig. 5 but the PCA is conducted using actual GOME-2 data taken over ice- and snow-covered surfaces, the Sahara, and 
cloudy ocean. 


optically thick and absorption is temperature dependent, sim- 
ulations are needed to evaluate how well our simplified ap- 
proach will work in this spectral region. Our approach is flex- 
ible in that different fitting windows may be utilized. As will 
be shown below, the O 2 A-band can in fact be completely re- 
moved from an extended fitting window without much if any 
loss of information content. 

To the extent that our approach is shown to be success- 
ful, there are several advantages. By using PCs, we do not 
need to estimate parameters affecting O 2 and H 2 O absorption 
such as surface pressure, temperature and water vapor pro- 
files, and cloud and aerosol parameters that affect nearly all 
satellite observations. In addition, we eliminate dependence 
upon a precise specification of the instrument response func- 
tion. The PCs may also capture instrument artifacts that are 
otherwise difficult to quantify. 

As in DOAS retrievals, we use the logarithm of sun- 
normalized radiance spectra for the PCA. We first compute 
PCs with the simulation training data set. For comparison, 
we also generate PCs using actual GOME-2 satellite radi- 
ance data. For the GOME-2 PCA, we use spectra from a sin- 
gle day (1 May 2007) consisting of observations over sea 
ice, snow/ice-covered land, the Sahara, and cloudy ocean for 
pixels with 0 o < 75°. For the cloudy ocean data, we compute 
the reflectance at 670 nm (p67o) and use observations only for 
P670 > 0.7. 


For both real and simulated data, we normalize the spec- 
tra with respect to a second-order polynomial fit to wave- 
lengths not significantly affected by atmospheric absorption 
(i.e., 712-713 and 748-757 nm, and wavelengths > 775 nm) 
before taking the logarithm of the spectra. This normaliza- 
tion produces values representative of the total sun to satel- 
lite absorptance. We are able to use a second-order polyno- 
mial for this fit of spectra that do not contain complex sur- 
face reflectances. However, as noted below, we need to use 
higher-order polynomials to characterize vegetated surfaces. 
Alternatively, PCAs may be similarly performed without tak- 
ing the logarithm of the normalized spectra in order to model 
transmittance instead of absorptance. 

Absorption affecting fluorescence in the far-red emission 
feature includes that from the O 2 A-band near 760 nm as well 
as a weaker water vapor band at shorter wavelengths. Fig- 
ure 4 shows examples of normalized spectra approximating 
the sun-to-satellite transmittance separately for two wave- 
length ranges: (1) 712-747 nm, dominated by H 2 O absorp- 
tion, and (2) 747-783 nm encompassing the O 2 A-band. 

Figure 5 shows the leading four PCs for the wave- 
length range 712-747 nm computed with simulated data for 
FWHMs of 0.5 nm (similar to GOME-2) and 0.3 nm. The 
spectral variance in this window is due almost exclusively 
to water vapor absorption. The variances explained (with 
respect to the total) as well as the cumulative variances 
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Fig. 7. Similar to Fig. 5 but for the long- wavelength fitting window using simulated data. 


explained are indicated. The PCs are similar for the two spec- 
tral resolutions with somewhat deeper structures at the higher 
resolution. The variance explained by the leading PCs is sim- 
ilar for the two spectral resolutions. 

Figure 6 similarly shows PCs generated from actual 
GOME-2 satellite data. PCs and variance explained are simi- 
lar for the simulated and GOME-2 data. The first PC explains 
over 99 % of the spectral variance, and 99.99 % of the vari- 
ance is captured in the first four modes for both the simulated 
and GOME-2 data. However, PC #4 from GOME-2 appears 
to correspond to PC #3 from the simulated data, but there 
is no similar correspondence between PC #3 from GOME-2 
and PC #4 from the simulation. The PCs for simulated and 
real data are expected to be different as PCs from the real data 
may contain information related to instrumental artifacts and 
processes not included in the simulated data (e.g., rotational 
Raman scattering). In addition, the simulated data may not 
represent all of the conditions or the distribution of condi- 
tions that are present in the GOME-2 data, particularly bright 
scenes that occur over heavily clouded conditions. The train- 
ing data set, however, does include bright soils and snow. 

Figures 7 and 8 similarly show the leading PCs for the 
spectral window 747-783 nm dominated by strong oxygen 
A-band absorption near 760 nm. Again, the PCs are similar 
for simulated and GOME-2 satellite data with over 99.9 % 
of variance captured by the leading mode, and more than 
99.999 % of the variance explained by the first four modes. 


PCs #2 and #3 appear to be reversed for the GOME-2 and 
simulated data. 

4.3 Solving the nonlinear problem 

To solve the nonlinear estimation problem, we use a gradient- 
expansion algorithm adapted from Marquardt (1963) and 
Bevington (1969). This algorithm provides a relatively fast 
convergence, typically 4-6 iterations. We derive and supply 
to this algorithm the analytic Jacobians or partial derivatives 
of the observed radiances with respect to the state variables. 

Typical Jacobians (i.e., partial derivatives of the re- 
flectances with respect to the coefficients of the PCs, sur- 
face reflectance polynomials, and the peak value of the far- 
red fluorescence feature at 736.8 nm) are shown in Fig. 10 
for FWHMs of 0.5 and 0.3 nm. Although the components 
are not completely orthogonal, our simulation results will 
show that fluorescence can be successfully disentangled from 
atmospheric and surface parameters. Subtle differences in 
the Jacobians enable this differentiation. For example, small 
Fraunhofer structures can be seen in the fluorescence Jaco- 
bian at wavelengths between about 745 and 758 nm that are 
not seen in the other Jacobians. An instrument with a high 
enough spectral resolution and SNR should be able to detect 
these features as will be demonstrated below. Again, deeper 
spectral structures are seen at the higher spectral resolution. 
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Fig. 8. Similar to Fig. 6 but for the long- wavelength fitting window computed using actual GOME-2 data. 
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Fig. 9. Flow diagram showing the overall processing scheme for GOME-2 fluorescence retrievals. 
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Fig. 10. Typical Jacobians {dp/dx) where v is the far-red fluorescence peak value (top), the wavelength-independent component of the 
surface reflectivity (middle), and the coefficient of the first PC (bottom) for FWHM = 0.5 nm (black) and FWHM = 0.3 nm (red). The PC 
analyses are carried out separately for the wavelength ranges 712-747 nm and 747-783 nm but are shown in a single plot for convenience. 


At convergence, the partial derivatives contained in the Ja- 
cobian K matrix may be used to compute errors from an un- 
constrained linear error estimation, i.e., 

S r =(K r S- 1 K) -1 , (7) 

where S r is the retrieval error covariance matrix, and S e is 
the measurement error covariance. We use this approach to 
compute fluorescence error standard deviations and correla- 
tions and compare these error estimates with those obtained 
in the full retrieval simulation as will be discussed below. We 
also use the linear approach to estimate errors from GOME-2 
retrievals. 

4.4 Processing of GOME-2 data 

The overall processing of GOME-2 data follows the flow di- 
agram shown in Fig. 9. One subset of radiance data is used 
to generate the PCs (chosen such that fluorescence is not 
present), and a different subset is used for the fluorescence 
retrievals (data over land). A PC A is conducted daily for each 
subwindow using the daily measured solar flux. We use data 
over highly cloudy ocean and snow- and ice-covered surfaces 
(p at 670 nm > 0.7) and the Sahara for the PCA. The derived 
PCs are then used for the fluorescence retrieval. Following 
the retrieval, quality assurance checks are done to remove 
highly cloudy data and failed retrievals as described in more 
detail below. Finally, the data are gridded at various tempo- 
ral resolutions to produce level 3 data sets. More than 5 yr of 



Fig. 11. Retrieved spectral surface reflectances for the range of con- 
ditions in the simulated data set. 


GOME-2 data have been processed, and the level 2 data are 
available from http://acdb-ext.gsfc.nasa.gov/People/Joiner/. 

4.4.1 GOME-2 fluorescence retrievals 

We use a fitting window from 712 to 783 nm for GOME-2 
retrievals with a fourth-order polynomial to model the sur- 
face reflectivity and a Gaussian function for the canopy- 
level far-red fluorescence as described above. We used 25 
PCs for each of the two PCA subwindows shown above. No 


Atmos. Meas. Tech., 6, 2803-2823, 2013 


www.atmos-meas-tech.net/6/2803/2013/ 


J. Joiner et al.: O 2 A-band fluorescence retrievals 

adjustments are made to the calibrated radiances/irradiances; 
it should be noted that the MetOp-A GOME-2 is known to 
have suffered from radiometric degradation over its lifetime. 

4.4.2 Cloud filtering and quality control 

As in Joiner et al. (2012), we compute the effective cloud 
fraction / c and eliminate data with f c > 0.4. To compute 
/ c , we use the black-sky 16-day gridded filled-land surface 
albedo product from Aqua MODIS (MOD43B3) at 656 nm 
(Lucht et al., 2000). Application of more or less stringent 
limits on cloud contamination within a moderate range (0- 
0.5) did not substantially alter the derived spatial and tem- 
poral patterns of F s , but stricter limits decrease the number 
of samples included, and therefore reduce coverage and in- 
crease noise in the gridded averages. 

In the results shown below, we include all data passing 
gross quality assurance checks on the retrieval convergence 
and radiance residuals. These checks removed few obser- 
vations in general. However, the South Atlantic Anomaly 
(SAA) adds noise to GOME-2 measurements in the vicin- 
ity of South America, and most of the data removed by the 
residual checks were located in this area. The distribution and 
sampling frequency of the data used to the generate the PCs 
leads to higher retrieval uncertainties than the SAA; the rel- 
atively small number of PCs used in the retrieval generally 
does not capture the highly variable errors found in the re- 
gion. While many of the affected spectra are screened out by 
radiance residual checks, the filtering process reduces sam- 
pling and may not eliminate all affected spectra. Therefore, 
GOME-2 errors are generally higher over South America as 
compared with other areas. We also eliminate all data with 
SZA > 70°. 

5 Sensitivity analysis 

In this section, we retrieve fluorescence using the simulated 
radiances for the 230400 different conditions contained in 
the testing data set described in Sect. 3. We conduct retrievals 
for a number of different scenarios. We then compare the re- 
trieved fluorescence with that of the truth as specified in the 
testing data set for the entire sample. Table 1 provides statis- 
tical results of those comparisons for the scenarios described 
below. 

5.1 Sensitivity to number of PCs used 

In the first experiment, we simulate data for a fitting win- 
dow between 747 and 780 nm and for an instrument with 
FWHM = 0.5 nm, sampling of 0.2 nm, and SNR = 2000. In- 
strument noise is uncorrelated between channels, and follows 
a Gaussian distribution. Here, we use a fourth-order polyno- 
mial to model the surface reflectivity. Increasing the polyno- 
mial order does not significantly improve the results, while 
use of a second-order polynomial significantly degrades re- 
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suits. For reference, Fig. 11 displays the range of spectral 
surface reflectances that are retrieved from the diverse simu- 
lated data. 

The first four lines of Table 1 show results for retrievals 
that use 5, 10, 15, and 25 PCs. There are small biases in all 
cases with biases generally decreasing with increasing num- 
bers of PCs. The improvement in both accuracy and precision 
is noticeable when increasing from 5 to 10 PCs, and levels 
out with further increases. There is virtually no change in the 
results when we increase the number of PCs from 25 to 35 
(not shown). We find that a similar number of PCs is needed 
for the shorter- wavelength subwindow (712-747 nm). For 
the remainder of this section, all retrievals will use 25 PCs 
for each of the two sub windows. 

We also computed fluorescence errors using the linear esti- 
mation method. Here, we assumed random and uncorrelated 
radiance errors. The radiance errors were also assumed to be 
independent of wavelength. This is consistent with the errors 
generated in the simulated data set. We computed errors for 
the various wavelength ranges, numbers of PCs, and SNRs 
used above. We find that the computed fluorescence error 
standard deviations from the linear approach are very close 
(within ~ 4 %) to those obtained in the full retrieval simula- 
tion. We find that fluorescence errors can be moderately cor- 
related with errors in the constant reflectivity term and some 
of the leading PCs. 

5.2 Sensitivity to signal-to-noise ratio 

For comparison, Table 1 shows results of fluorescence 
retrievals using simulated observations as above with 
SNR = 1000. Results may be compared with the 25 PC 
case (line 4) described in Sect. 5.1 for an instrument with 
SNR = 2000. Standard deviations and root-mean- squared 
(RMS) errors for SNR= 1000 are slightly less than twice 
those at SNR = 2000. This is consistent with the expected 
behavior of a retrieval based solely on solar Fraunhofer line 
structure, where errors are inversely proportional to SNR 
with all other parameters held constant (Joiner et al., 2012). 

5.3 Sensitivity to the fitting window 

In the experiments shown in lines 6-10 of Table 1, we use 
different fitting windows for an instrument with SNR = 2000 
and FWHM = 0.5 nm. Fluorescence errors are approximately 
a factor of 2 smaller with the 712-747 nm fitting window as 
compared with the 747-780 nm shown above. This demon- 
strates that fluorescence retrievals can be obtained with good 
precision without using the filling-in signal from the O 2 A- 
band. The 712-747 nm window spans the peak far-red fluo- 
rescence, and also contains H 2 O absorption. A reduced fit- 
ting window of 723-747 nm produces slightly degraded re- 
sults as compared with the 71 2-747 nm window. Both of 
these shorter-wavelength windows show decreased bias as 
compared with the 747-780 nm window that includes the O 2 
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Table 1. Statistical comparison of retrieved versus true values of F s obtained with the simulated testing data set for different experiments 
(Exp); all fluorescence radiance units (indicated by *) are mWm -2 nm -1 sr _1 . Retrievals are performed for an instrument with a given 
full-width-at-half-maximum (FWHM) line shape function, signal-to-noise ratio (SNR), number of principal components (#PCs), and fitting 
window from starting wavelength to ending wavelength X 2 . Statistics given are the root-mean- squared difference (RMS diff.), correlation 
coefficient (r), mean difference (bias) of retrieved minus truth, standard deviation (<r), and slope (A) and intercept (B) of a linear fit (retrieved 
fluorescence = A + B • truth) (last two columns). 


Exp 

FWHM 

(nm) 

SNR 

#PCs 

V 

(nm) 

*2 

(nm) 

RMS diff. 

* 

r 

bias 

* 

0 

* 

slope 

intercept 

* 

1 

0.5 

2000 

5 

747 

780 

0.88 

0.67 

-0.68 

0.57 

0.43 

-0.03 

2 

0.5 

2000 

10 

747 

780 

0.58 

0.82 

-0.37 

0.45 

0.58 

0.11 

3 

0.5 

2000 

15 

747 

780 

0.50 

0.87 

-0.32 

0.38 

0.72 

-0.01 

4 

0.5 

2000 

25 

747 

780 

0.43 

0.87 

-0.22 

0.38 

0.80 

0.01 

5 

0.5 

1000 

25 

747 

780 

0.70 

0.69 

-0.22 

0.67 

0.80 

0.01 

6 

0.5 

2000 

25 

712 

747 

0.33 

0.97 

-0.07 

0.32 

1.00 

0.07 

7 

0.5 

2000 

25 

723 

747 

0.37 

0.97 

-0.01 

0.37 

0.94 

0.13 

8 

0.5 

2000 

25 

755 

775 

0.56 

0.76 

-0.30 

0.47 

0.71 

0.01 

9 

0.5 

2000 

25 

759 

768 

1.19 

0.43 

-0.46 

1.10 

0.70 

-0.13 

10 

0.5 

2000 

25 

747 

758 

1.48 

0.58 

-0.03 

1.48 

0.98 

0.01 

11 

0.3 

2000 

25 

747 

780 

0.49 

0.85 

-0.29 

0.40 

0.80 

-0.06 

12 

0.3 

2000 

25 

712 

747 

0.22 

0.99 

-0.11 

0.19 

0.94 

0.02 




Fig. 12. Fluorescence retrievals from simulated data (y axis) using wavelengths between 712 and 747 nm for instruments with 
FWHM = 0.5 nm (left) and FWHM = 0.3 nm (right), both with a signal-to-noise ratio of 2000. Fluorescence is averaged over the wave- 
lengths used in the retrieval and compared with the “truth” (x axis) averaged in the same way. Standard deviations are shown with vertical 
bars. Different symbols are shown for the various values of chlorophyll content, and different colors are for the different values of leaf area 
index. 


A-band. Our results are consistent with those of Guanter et al. 
(2013), who similarly showed that fitting windows in this 
spectral region, without the benefit of the O 2 A-band, can be 
used to retrieve fluorescence with higher-spectral-resolution 
ground-based instruments. 

Table 1 also examines results for fitting windows more 
confined to the O 2 A-band. As may be expected, retrievals 
are degraded for a smaller fitting window of 755-775 nm as 
compared with the 747-780 nm window used above. Results 
are also shown with a smaller fitting window of 759-768 nm. 
This window consists primarily of the O 2 A-band with em- 
bedded Fraunhofer structure. Note that the strongest solar 
Fraunhofer feature within the range 712-783 nm occurs in- 
side the O 2 A-band region near 766 nm. Results are signif- 
icantly degraded with this limited fitting window. Franken- 


berg et al. (201 la) noted the difficulties associated with dis- 
entangling fluorescence spectral information from that of 
aerosols, clouds, surface pressure, etc., using only wave- 
lengths within O 2 A-band spectral region. Although there is 
a significant correlation between fluorescence and these other 
parameters, there is nevertheless a limited ability to retrieve 
information about fluorescence within this absorption band. 
However, biases increase as the fitting window is more con- 
fined to the O 2 A-band spectral region. Precision is signifi- 
cantly degraded with a limited fitting window of 747-758 nm 
containing only weak solar Fraunhofer line structures, while 
the accuracy for this fitting window is good. 
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Wavelength (nm) 

Fig. 13. RMS of simulated radiance residuals (in percent of 
radiance) from the testing data set with FWHM = 0.3 nm and 
SNR = 2000 when fluorescence radiance (F s ) is fitted/retrieved and 
when it is not fitted/retrieved. 


5.4 Sensitivity to spectral resolution 

In order to show the sensitivity of the results to spectral 
resolution, we performed a similar set of experiments at 
a higher spectral resolution (FWHM = 0.3 nm, sampling of 
0.1 nm) as compared with that used in the above subsec- 
tions. Lines 12 and 13 of Table 1 show retrieval statistics for 
the higher-spectral-resolution instrument with SNR of 2000 
for fitting windows of 747-780 nm and 712-747 nm, respec- 
tively. The precision is significantly improved as compared 
with FWHM = 0.5 nm retrievals. This improvement results 
from (1) more spectral samples within the fitting window and 
(2) a larger filling-in fluorescence signal in the cores of the 
deeper Fraunhofer lines that are better resolved at the higher 
spectral resolution as shown in Fig. 10. Further improve- 
ments can be made by making measurements with higher 
spectral resolution and/or increased sampling. We also var- 
ied the fitting window as in Sect. 5.3 for FWHM = 0.3 nm, 
and reached the same conclusions as with FWHM = 0.5 nm. 

As described above, there are 60 distinct values of fluo- 
rescence (averaged between 740 and 780 nm) ranging from 
near zero to near 4mWm -2 nm -1 sr -1 . For each value, the 
observing conditions vary (e.g., different SZAs, VZAs, sur- 
face pressures, temperature profiles, and aerosol parameters). 
Figure 12 shows retrieval results using the 71 2-747 nm fit- 
ting window for the FWHM = 0.3 and 0.5 nm simulated data 
with SNR = 2000. This figure shows that biases are more 
prevalent for higher levels of fluorescence. 

5.5 Radiance residuals from simulated data 

Figure 13 shows RMS of the radiance residuals (ob- 
served minus calculated radiance) for two cases with 
FWHM = 0.3 nm, SNR = 2000, and a fitting window of 747- 
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Fig. 14. RMS of GOME-2 radiance residuals obtained with and 
without fitting F s (in percent of the observed radiance) for a sin- 
gle day for moderately to highly vegetated pixels. 

780 nm: (1) fluorescence is not retrieved, and (2) fluores- 
cence is retrieved. The RMS of the residual is computed at 
each wavelength averaged over all conditions in the simula- 
tion testing data and shown as a percentage of observed radi- 
ance. As expected, reductions in the residuals are achieved 
when fluorescence is retrieved, particularly in the vicin- 
ity of deep solar Fraunhofer lines (e.g., near 749, 751, 
and 766 nm). These improvements occur throughout the full 
spectral range. 

We conducted an additional experiment in which we use a 
flat solar spectrum in place of the simulated solar spectrum 
with the fitting range 747-780 nm that includes the O 2 A- 
band. We obtained poor results with our retrieval approach 
when using the flat solar spectrum (essentially no informa- 
tion about F s is retrieved). This is another indicator that our 
method relies primarily on the filling-in of solar Fraunhofer 
lines to retrieve F s . 

Note that relatively larger residuals (larger than instrument 
noise but well below 1 %) are produced at the very low radi- 
ance levels found within the deep O 2 A-band. This result is 
seen in both the simulated data as well as real GOME-2 data 
as will be shown below. These residuals could be the result 
of nonlinear behavior of the O 2 A-band that the PCA method 
is not able to capture. 

5.6 Discussion of simulation results 

As noted above, our retrieval approach relies on several sim- 
plifying assumptions. For example, we assume that atmo- 
spheric scattering is negligible, that the radiative transfer as- 
sumptions are only valid for monochromatic light, and that 
the spectral structures of F s and p s could be modeled with 
a few parameters. The simulated data contain none of these 
assumptions; the radiances are generated monochromatically 
with scattering before they are convolved with the instrument 
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(a) GOME-2 (left) and GOSAT (right), July 2009 



(b) GOME-2 (left) and GOSAT (right), December 2009 



(c) GOME-2 (left) and GOSAT (right), Annual 2009 


Fig. 15. Global composites of F s from GOME-2 (left) and GOSAT-FTS (right) retrievals for July (top), December (middle) and the annual 
(bottom) average in 2009 (June 2009 through May 2010 for GOSAT). GOME-2 are binned in 0.5° cell boxes; GOSAT retrievals binned in 
2° cell boxes. Both retrievals refer to wavelengths near 757 nm. 


response function, and the spectral dependences of F s and 
p s are based on model and spectral libraries. Therefore, our 
simulation results should accurately reflect errors produced 
by these assumptions. As can be seen, the biases and errors 
produced by these simplifications are relatively small. We did 
not simulate RRS in the simulation. Our assumption in pro- 
cessing GOME-2 data is that the PCA will be able to disen- 
tangle the spectral effects of RRS from those of F s . This will 
be discussed further below. 

6 Results from GOME-2 data 

6.1 Radiance residuals from GOME-2 

Figure 14 shows the spectral RMS of the radiance residu- 
als, obtained with and without fitting F s , over the GOME- 
2 fitting window for a single day. Residuals are averaged 


for each wavelength over all observations with SZA <70° 
and NDVI > 0.3 (i.e., moderately to highly vegetated pixels) 
that passed quality-control and cloud-filtering checks. Re- 
ductions in the residuals can be seen throughout the spectral 
range with a similar spectral structure as shown in Fig. 13 for 
the simulated data (i.e., reductions at the deep solar Fraun- 
hofer lines). Some reduction is also seen at the deepest part 
of the O 2 A-band near 760 nm. In the portions of the fitting 
window that are relatively free of atmospheric absorption, 
the residuals are consistent with a GOME-2 SNR of 1 000 or 
greater. Residuals have similar magnitudes to those shown in 
Fig. 13 for simulated data with SNR = 2000. Residuals are 
somewhat higher in regions where atmospheric absorption is 
present such as in the water vapor absorption band shortward 
of about 747 nm. 
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Fig. 16. Global maps of GOME-2 F s retrieval statistical parameters in a 0.5° grid cell for July (left column) and December (right column) 
2009. Each column shows the standard deviation (top) and the number of points per grid cell (bottom). 


6.2 Comparison of GOME-2 and GOSAT fluorescence 

Global composites of F s derived from GOME-2 for July, 
December and the annual average in 2009 are displayed in 
Fig. 15. Here, we converted the retrieved peak values of flu- 
orescence at 736.8 nm to values that would have been ob- 
served at the GOSAT wavelengths (near 759 nm). We multi- 
plied the 736.8 nm fluorescence by a factor of 0.59. This fac- 
tor is consistent with the spectral shape of the fluorescence 
emission assumed in the retrievals. It is also consistent with 
the fluorescence emissions used in the simulated testing data 
set (see Fig. 3). For comparison, maps of F s at 757 nm from 
GOS AT-FTS retrievals are also shown for the same time peri- 
ods (the annual average is from June 2009 through May 2010 
for GOSAT). The GOSAT-FTS retrievals were performed 
and processed with the algorithm described in Guanter et al. 
(2012). Quality-filtered GOME-2 retrievals have been aver- 
aged in 0.5° latitude-longitude grid boxes, whereas a 2° grid 
is used for the GOSAT retrievals that have a much more 
sparse spatial sampling than GOME-2. The filtering uses data 
only for SZA< 70°, where RRS effects should be small. We 
see no obvious biases resulting from RRS for these condi- 
tions. 

A very good agreement of the F s spatial patterns is ob- 
served between the two data sets, although improvements in 
spatial resolution and precision are obvious in the GOME-2 
maps. High F s values are observed over densely vegetated 
areas. Globally, the highest F s signal is found in July in 
the eastern United States. High F s values are also observed 


in some parts of South America and Africa in December. 
In some of these regions, GOME-2 values are higher than 
GOSAT. This may be the result of differences in absolute 
calibration of the instruments or small biases in the GOME-2 
algorithm. The higher GOME-2 values are inconsistent with 
the differences due to time of day and cloud contamination 
as discussed in more detail below. 

F s values near zero are detected in July over Greenland 
(zero to slightly negative), in December over Antarctica (zero 
to slightly positive), and during the entire year over the 
Sahara and most of Australia. These spatial patterns com- 
pare well with those observed in the F s maps derived from 
GOSAT-FTS for the same time periods. The slight biases 
(both positive and negative) in these areas are likely not re- 
lated to random instrumental noise as this would be removed 
in a long-term average. There are several potential sources of 
these small biases, including simplifying assumptions in the 
forward model, small correlations between the fluorescence 
spectral signal and that of reflectivity and/or some of the PCs, 
lack of representativeness of the PCs used in the retrieval, and 
systematic instrumental artifacts. 

Concerning the annual averages of F s retrievals, the main 
difference between the two data sets is in the tropical rain- 
forest, especially in Africa and Indonesia. GOME-2 F s is, in 
relative terms, lower than that of GOSAT over those areas. 
This may be due to a larger impact of cloud contamination in 
coarser GOME-2 footprint data (40 km x 80 km for GOME- 
2 as compared with an around 10 km diameter for GOSAT). 
The smaller footprint of GOSAT may also allow for higher 
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Fig. 17. Ten-day composites of F s and NDVl derived from GOME-2 data in 0.5 x 0.5 0 grid cells between day of year (DOY) 131 and 160 
of 2009. Fluorescence is normalized by the cosine of the sun zenith angle (/x s ) in order to minimize the temporal and latitudinal dependence 
of fluorescence on incoming at-surface photo synthetically active radiation. 


peak values of F s signals to be obtained in areas of very dense 
productive vegetation. At the same time, the low sampling of 
GOSAT may not allow for adequate representation of grid- 
box averages. 

The different solar illumination angles encountered by the 
two instruments may also contribute to relative differences 
between high and low latitudes. The overpass time of the 
satellites (~ 09:30 LT for MetOp-A, - 13:00 LT for GOSAT) 
results in GOME-2 measurements that have systematically 
higher solar zenith angles (SZAs) as compared with GOSAT. 
In general, the illumination angle affects the fluorescence sig- 
nal at the top of canopy through (1) the intensity of the sun- 
light incident at the canopy; (2) the amount of illuminated 
leaves that is related to the ratio of of diffuse-to-direct irradi- 
ance; and (3) the physiological relationship between photo- 
synthesis, fluorescence, and heat dissipation. 


The overall excellent spatial agreement between GOME-2 
and GOSAT can be considered as a rough first validation of 
our GOME-2 F s retrievals given the fact that the GOME-2 
retrieval approach is much more complicated and prone to 
systematic errors as compared with that of GOSAT. While 
the GOME-2/GOSAT comparison may be considered as a 
consistency check and the simulation experiments provide 
evidence that the retrieval approach is valid, there are short- 
comings to these types of evaluations. Firstly, the GOSAT 
retrievals themselves have only been evaluated using indi- 
rect strategies such as plausibility checks and simulation ex- 
periments. Secondly, the simulation studies, while quite de- 
tailed, generally do not contain all the complexities found in 
real satellite data (e.g., instrumental artifacts and RRS). A 
more direct validation using ground- and aircraft-based data 
remains challenging owing to the large pixel sizes of the cur- 
rent satellite instruments. Approaches need to be developed 
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and tested to scale up measurements made at small scales 
(ground-, aircraft-, and small-footprint satellite data) to the 
larger GOME-2 pixels. 

6.3 Estimated errors, variability, and numbers of 
GOME-2 retrievals 

The standard deviation of the July and December F s re- 
trievals in Fig. 16a indicates the variability of the F s values 
observed in each 0.5° x 0.5° gridbox. Sources of this variabil- 
ity are instrumental noise, natural variability in vegetation 
activity within the month, residual cloud effects, variabil- 
ity owing to different illumination and viewing geometries, 
and the different footprints from various MetOp-A orbits. 
We have observed systematic variability with respect to view 
zenith angle for a given area over the course of a month as 
has been reported previously for GOSAT and SCIAMACHY 
data (Guanter et al., 2012; Joiner et al., 2012). The analysis 
of the effects of the illumination and viewing geometries on 
the F s signal will be explored in future works. 

It can be seen that the highest variability in the retrievals is 
found over a large area in South America, and that this does 
not depend on the season. Instrument performance in this re- 
gion is substantially degraded by the SAA. Even though this 
effect does not appear to have a large impact on the monthly 
average in Fig. 15a and b, the data over this area must be han- 
dled carefully. Concerning the rest of the globe, the standard 
deviation patterns compare well with the expected at-sensor 
radiance patterns (e.g., higher standard deviations over bright 
snow-covered areas and deserts) that can be explained by the 
higher contribution of photon noise triggered by higher at- 
sensor radiance levels. 

We estimated errors based on the linear approach de- 
scribed above assuming random and uncorrelated radi- 
ance errors with a signal -to-noise ratio of 1000 consistent 
with the radiance residual analysis. We obtained fluores- 
cence standard deviation values similar to the variability 
shown in Fig. 16 (estimated standard deviations of ~0.4- 
1 mW m -2 nm -1 sr -1 in moderately to heavily vegetated ar- 
eas). However, the linear error estimates were higher than 
the variability shown over surfaces with high reflectivity in- 
cluding Greenland and the Sahara region. Estimated errors 
in these areas can exceed 1.5 mW m -2 nm -1 sr -1 . This may 
indicate that the SNR of GOME-2 is higher for these bright 
scenes as compared with other areas. Error estimates were 
lower than the variability shown in the region of the SAA, 
because the increased noise in this area was not accounted 
for in the linear estimation. 

Figure 16b shows the number of F s retrievals fulfilling the 
quality criteria per gridbox. Note that the red orbital stripes 
are not artifacts but are due to narrow swath data that are 
obtained approximately once per month. Typically, 10-25 
retrievals are available for each grid box within a month. 
Those numbers are smaller over highly cloudy tropical rain- 
forest regions, especially during the wet season in Decem- 
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ber, and at high latitudes. Therefore, fluorescence standard 
errors for the monthly mean gridded data are reduced by fac- 
tors between approximately 3 and 5 as compared with the 
estimated single observation error. Estimated fluorescence 
errors for monthly mean gridded data are in the range 0.1- 
0.4 mW m -2 nm - 1 sr - 1 . 

6.4 Temporal variations in GOME-2 fluorescence and 
NDVI 

The high revisit time of GOME-2 allows for excellent tem- 
poral sampling in the derived vegetation products. This is 
illustrated in Fig. 17. The figure shows 10-day composites 
of F s and NDVI derived from GOME-2 data between day 
of year (DOY) 131 (11 May) and 160 (10 June) of 2009. 
Here, fluorescence values are normalized by the cosine of 
the solar zenith angle in order to minimize the latitudinal and 
temporal variations in fluorescence owing to the incoming 
PAR. A lower signal-to-noise ratio is observed for these 1 0- 
day composites as compared with the monthly averages in 
Fig. 15, especially in the area of the South American region 
affected by the SAA. Spatial gaps in the data are due to per- 
sistent cloud contamination, as MetOp provides near-daily 
global coverage. 

We use a standard definition to compute NDVI from 
GOME-2, i.e., 

NDVI = PN ' R ~ PRED , (8) 

Pnir + Pred 

where the pnir and pred are computed using single wave- 
length observations closest to 780 and 670 nm, respectively. 
Note that the values of pnir and pred have not been cor- 
rected for atmospheric scattering, surface BRDF effects, or 
fluorescence, and are affected by cloud contamination within 
the GOME-2 footprint. Despite the simplicity of the GOME- 
2 NDVI calculation, spatial patterns are similar to those of 
the MODIS NDVI product (not shown). The GOME-2 NDVI 
sampling is identical to that of the GOME-2 fluorescence. 

Phenological changes in the Northern Hemisphere are 
clearly visible from one 10-day period to another in Fig. 17. 
A strong increase in F s is observed in Europe from DOYs 
131-140 to DOYs 141-150. This rapid change in green 
biomass is also detectable in the NDVI, although with 
a smaller intensity. This high temporal sampling of F s trends 
cannot be achieved with GOSAT data owing to the signifi- 
cantly lower number of observations and the sparse spatial 
sampling. Improved temporal sampling should be achievable 
by processing data from both GOME-2 instruments for the 
periods of dual measurements. 

7 Conclusions 

We have developed a new approach to retrieve far-red flu- 
orescence from moderate-spectral-resolution satellite instru- 
ments. The method utilizes fluorescence filling-in of the O 2 
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A and water vapor bands as well as the surrounding weak so- 
lar Fraunhofer lines; it relies upon the separation of spectral 
signatures produced by upwelling chlorophyll fluorescence 
and atmospheric absorption as well as surface, cloud, and 
aerosol backscattering of solar radiation. We use principal 
components, derived from data free of fluorescence, to esti- 
mate the spectral structure of atmospheric absorption. This 
information is incorporated into a simplified radiative trans- 
fer model that accounts for atmospheric absorption of fluo- 
rescence emissions. Through detailed simulations, we show 
that high-quality fluorescence retrievals can be obtained us- 
ing instrumentation with high SNR and moderate spectral 
resolution similar to GOME-2. Retrieval errors depend upon 
the instrument SNR, spectral resolution, and specification of 
the spectral fitting window. 

We then applied our new approach to satellite moderate- 
spectral-resolution measurements from GOME-2. The 
GOME-2 retrievals compare well with those from GOSAT 
that are processed with a less complex algorithm, providing 
further confidence in our approach and implementation with 
real data. Fluorescence errors for monthly mean gridded data 
are estimated to be ~ 0. 1-0.4 mW m 2 nm 1 sr _1 . Owing to 
the excellent spatial sampling and high signal-to-noise ra- 
tio of the GOME-2 measurements, we are able to map far- 
red terrestrial fluorescence at higher spatio-temporal resolu- 
tions than previously published GOSAT and SCIAMACHY 
data. This mapping shows clearly for the first time a north- 
ward shift in PAR-normalized fluorescence within the single 
month of May as the sun shifts northward during the boreal 
spring. 

Several satellite instruments with NIR spectral coverage 
and various spectral and spatial resolutions have flown, are 
currently flying, or are planned for launch in the next few 
years. The approach outlined here can potentially be ap- 
plied to these instruments. SCIAMACHY provides observa- 
tions in the same spectral region. While the native SCIA- 
MACHY footprint (30 km x 60 km) is slightly smaller than 
that of GOME-2 (40 km x 80 km), the spatial sampling of 
SCIAMACHY is not as good as GOME-2, in part due to 
alternating between limb and nadir measurements. In addi- 
tion, SCIAMACHY observations in the near-infrared at some 
wavelengths were spatially coadded and are not provided 
at full spatial resolution in the level lb data set. The orig- 
inal GOME instrument, launched in 1995 on the European 
Space Agency’s European Remote Sensing satellite 2 (ERS- 
2), can also be used for fluorescence measurements, but with 
a larger pixel size (40 km x 320 km) in its nominal operat- 
ing mode. SCIAMACHY and GOME have the unique abil- 
ity to extend the record of fluorescence measurements back 
to 1995. In addition to GOME-2 and SCIAMACHY, the ap- 
proach may also be applied to the GOSAT interferometer, 
the Orbiting Carbon Observatory-2 (OCO-2) (Crisp et al., 
2004), planned for launch in 2014, and the TROPOspheric 
Monitoring Instrument (TROPOMI) (Veefkind et al., 2012) 
to be launched in 2015. The FLuorescence EXplorer (FLEX) 
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(Rascher, 2007; European Space Agency, 2008), an ESA Ex- 
plorer 8 mission, selected for Phase A/B 1 in early 2011, plans 
to utilize the O 2 A- and B-bands for chlorophyll fluorescence 
retrievals (Guanter et al., 2010) and other bio-spectral infor- 
mation across the visible through NIR spectral range. FLEX 
would provide measurements at a higher spatial resolution 
than current satellite sensors that were not designed for fluo- 
rescence measurements. 
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